Bose-Einstein Condensation in a Trap: the Case of a Dense Condensate 



Klaus Ziegler and Alok Shukla 
Max-Planck-Institut fur Physik Komplexer Systeme, Auflenstelle Stuttgart, Postfach 800665, D-70506 Stuttgart, Germany 

(February 7, 2008) 



o 



G 

O 

o 



> 

oo 
m 
o 

o 

On 
-i— > 



i 

c 

o 
o 



X 



We consider the Bose-Einstein condensation of atoms in a 
trap where the density of particles is so high that the low den- 
sity approach of Gross and Pitaevskii will not be applicable. 
For this purpose we use the slave boson representation which 
is valid for hard-core bosons at any density. This description 
leads to the same results as the Gross-Pitaevskii approach in 
the low density limit, but for higher densities, it predicts the 
depletion of the condensate in the regions where the density 
of the atomic cloud is high. 
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I. INTRODUCTION 

Atoms in a magnetic trap present an interesting sys- 
tem for the analysis of highly degenerate gases. The trap 
plays the role of a three-dimensional confining potential 
well and by using advanced experimental techniques such 
as laser and evaporative cooling, it is possible to study 
the the gas over a wide range of parameters like tem- 
perature or density of particles. One of the most spec- 
tacular achievements of such techniques was the obser- 
vation of the Bose-Einstein (BE) condensation in gases 
composed of alkali atoms M. BE condensation, origi- 
nally studied in terms of the ideal (non-interacting) Bose 
gas requires a minimal density of bosons po = <?3/2(l)/A 3 , 
where A = (2-kTi 2 /inksT) 1 / 2 is the de Broglie wave length 
and 33/2(1) sa 2.612. On the other hand, it is known that 
the density of the condensate in an interacting Bose gas 
at high density is depleted if the total density exceeds a 
certain value. For instance, the condensate in the bulk 
of 4 He is only about 10% at zero temperature. On the 
surface, however, the condensate can reach almost 100% 
because of the reduced total density. This phenomenon 
was observed in numerical simulations of an interacting 
Bose gas Q and in analytic calculations including an at- 
tractive interaction 2] or in a slave boson approach to a 
hard-core Bose gas The effect can be understood as 
a reduction of long-range correlations, necessary for the 
formation of a condensate, which is caused by increasing 
fluctuations due to an increasing density of interacting 
particles. There is an approach to the dilute interacting 
Bose gas due to Ginzburg and Pitaevskii || and Gross 
H , analogous to the Ginzburg-Landau approach for sec- 
ond order phase transitions. As in the general Ginzburg- 
Landau approach, the Gross-Pitaevskii (GP) approach is 
an expansion in powers of the order parameter field up 
to fourth order. It works very well close to the critical 



point, however, away from it, where the order parameter 
is not small anymore, it may significantly deviate from 
the correct result. This is not a problem in a homoge- 
neous system with a uniform order parameter, since we 
can restrict the theory to a regime where the order pa- 
rameter is small. In an inhomogeneous system, e.g., in 
a trap, the order parameter varies in the system. There- 
fore, it is not sufficient, unless the system is very dilute, 
to assume that the order parameter is small in some spa- 
tial region because it can be large in another region of 
the system. The interparticle interactions in the GP ap- 
proach are approximated by a hard core two-body po- 
tential with contribution only from the s wave scattering 
length. This approximation is quite satisfactory for the 
extremely dilute gases satisfying the condition na 3 <C 1, 
where n is density of particles and a is the scattering 
length. The physical implication of this condition is that 
it is highly improbable for three or more particles to in- 
teract with each other simultaneously. Therefore, it is 
clear that as the density of particles in the condensate 
increases, the likelihood of the three-body and higher or- 
der interactions will also increase, making it necessary 
to go beyond the GP approach. Even though the recent 
experiments on the BE in magnetic traps were based on 
dilute systems of bosons with na 3 as 10~ 6 , it is foresee- 
able that experiments can be performed where the Bose 
gas is dense. This is already indicated by the history of 
these experiments: the number of particles in the con- 
densate has increased by three orders of magnitude Q| 
as compared to the early results. In this context it is in- 
teresting to study the condensate including higher order 
interaction effects, expected in systems at higher density. 
The purpose of this article is to apply a method, which 
goes beyond the GP approach, to analyse an interacting 
Bose gas in a trap at arbitrary densities. 

The rest of the paper consists of two parts: in the first 
(Sect. II) part we discuss the case of a dilute Bose gas 
using the GP approach. This includes a mean-field the- 
ory based on the non-linear Schrodinger equation and the 
Thomas-Fermi approximation. In the second part (Sect. 
Ill) the Bose gas with hard-core interaction is defined as a 
functional integral in a slave boson representation. From 
the latter we derive an effective functional integral for the 
order parameter field which describes the BE condensate. 
The new effective functional integral, which also consti- 
tutes the main result of this work, is valid for an arbitrary 
condensate density, and takes account of the three-body 
and higher order effects in the interparticle interactions, 
at a finite temperature. We apply this approach to the 
problem of BE condensation in a trap in Sec III A, and 
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study the behavior of the condensate wavefunction as a 
function of chemical potential (which controls the par- 
ticle density) by means of a mean-field theory based on 
the Thomas-Fermi approximation again. However, like 
the GP approach, the limitation of the present work is 
its inability to account for the effects of the atoms out- 
side of the condensate on the condensate itself. Such an 
extension will be the subject of a future publication. 



II. DILUTE BOSE GAS 

The Bose gas, defined as a grand canonical ensemble of 
bosons, can be described in second quantization, for in- 
stance, using a functional integral representation H . The 
fluctuations with respect to time (i.e., quantum fluctua- 
tions) have a gap fc^T at non-zero temperature T due 
to the Matsubara frequencies u„ = rifc^T (n = 0, 1, ...). 
The condensation of bosons is characterized by a spon- 
taneous breaking of a [/(l)-symmetry (the phase degree 
of freedom of the complex boson field). This implies a 
Goldstone mode which describes gapless fluctuations in 
space. Thus the latter fluctuations are relevant for the 
condensation whereas the fluctuations with w„ ^ can 
be neglected because of the gap ksT. Thus it is suffi- 
cient to consider the wo-component of the quantum field 
$(x, wo) = since we are only interested in static prop- 
erties of the condensate near the phase transition. This 
approximation has been used in the discussion of the con- 
densation of a Bose gas in translational invariant systems 
[^],[| as well as in a harmonic trap P-fT 7 ^ . 

As an introduction we present the GP approach of a 
grand canonical Bose gas. The latter is defined by the 
partition function M 



Z = J e- Sop Y[d$ x d$* x . 



(1) 



The bosons in the condensate are described by the com- 
plex field $ which is controlled by the GP action of a 
dilute hard-core Bose gas 

Sgp = £ ■!',/, ,■<!-;■ - ]>>l^| 2 - ^l 4 ). (2) 



Since this model is based on a Bose gas with hard-core 
interaction [fl3|| , only three independent parameters en- 
ter: the chemical potential fi, the scattering length of 
the hard-core interaction a and the mass of the bosons. 
The coupling constant u is proportional to the scatter- 
ing length p3| , whereas the mass enters into the hopping 
rate r = ^ t x>x i < which has the same energy scale 
as /i. Formally, the action (2) is defined on a lattice with 
lattice constant a. This is a good approximation of the 
hard-core gas with scattering length a if one is only in- 
terested in length scales large compared to a. This is the 
case in the experiments because the typical size of the 
condensate is 5 x 10~ 4 cm whereas a«5x 10 -7 cm [H. 



Usually the hopping term J2 X x > $xt x ,x'$ x ' is replaced 
by the continuum approximation $ x r(l + (l/6a 2 )V 2 )<I>* 
and the sum by a formal integral, for simplicity. 

The magnetic trap can be modeled by introducing a 
confining potential V x , for instance, a harmonic potential. 
It is convenient to use a dimensionless expression for the 
action in (2) as given in Refs. |) 10|. The parameters 
of a gas with about 5000 87 Rb atoms, studied in the 
experiment by Anderson et al. Q, has been estimated in 
Ref. jl(J: The shape of the trap is anisotropic with the 
potential V x — x\ + x\ + 8ir§ — (j,q. The coupling constant 
of the interaction of the 87 Rb-atoms is u « 813 and the 
effective chemical potential /iq — t ~ 16.3. 



A. Mean-Field Theory 

The properties of the condensate (e.g., the density) 
can be evaluated using the saddle point approximation 
of the action Sgp- This is equivalent to the approxi- 
mation which neglects fluctuations of & x (classical field 
approximation). The classical field $ x is a solution of 
the non-linear Schrodinger equation 



((r/6a 2 )V 2 



t + V x + u|$ x | 2 )$ x = 0. 



(3) 



This nonlinear differential equation is complicated and 
generally one has to resort to numerical methods to ob- 
tain its exact solutions. However, in order to understand 
its behavior in the high-density limit one can neglect the 
kinetic energy term because, in that case, the nonlinear 
term of the equation is dominant. This is known as the 
Thomas-Fermi approximation. In a translational invari- 
ant system (i.e., V x = —po) the Thomas-Fermi approxi- 
mation gives a linear behavior for the condensate density 
as a function of fj,Q 



Pc oc l^l 2 = (1/m)(mo - t)6(^ - t), 



(4) 



where is the Heaviside step function. The condensate 
density increases ad infinitum upon increasing fiQ. This 
behavior is in disagreement with the depletion of the con- 
densate expected at higher total densities. It also con- 
tradicts the fact that the hard-core potential limits the 
density of the Bose gas. The reason behind this behav- 
ior is that in the GP approach the hard-core condition 
is implemented by a two-body delta function potential. 
However, in the high-density limit when a large number 
of particles are close to each other, this potential does not 
provide a strong enough repulsion needed for the strict 
imposition of the hard-core condition. In other word, cor- 
rections due to three-body and higher order effects will 
become equally important. This reflects that the GP ap- 
proach is realistic only for a dilute Bose system: The 
linear behavior is correct near the critical point where 
$| 2 ~ but the slope of the density of the condensate is 
less than linear as one goes away from the critical point. 
And finally it decreases, indicating the depletion of the 



2 



condensate at higher total densities. By the same argu- 
ment one can also obtain the inhomogeneous condensate 
density for the case of a high-density Bose gas composed 
of alkali atoms in a trap confined by a potential V x H 



Px cx l^l 2 = -(l/u)(V x + t)Q(-V x - r). 



(5) 



One may try to correct the high-density behavior of GP 
approach by including terms such as l^l 6 , I^zl 8 , - • ■ in 
the GP action (Eq.(g)) to account for the three-body and 
higher order interactions terms. But our feeling is that 
one will need to go to very high orders in this expansion 
to obtain the correct limiting behavior. However, in this 
work we propose an alternative approach which restores 
the correct high-density density behavior of an inhomo- 
geneous Bose gas by imposing strict hard-core condition 
by adopting a slave boson representation. This approach 
is an extension of a similar approach formulated for the 
case of a homogeneous Bose gas earlier, by one of us [||. 
In the following section we will develop and apply the 
aforementioned slave boson based approach to study the 
hard-core Bose gas in a harmonic trap. 



III. SLAVE BOSON APPROACH 

The slave boson representation was originally devel- 
oped for fermion systems (e.g., Hubbard model) The 
advantage of using a slave boson representation in case 
of strongly interacting electrons is that one can account 
for many effects of strong correlations at the mean-field 
level of this representation Q. For the Bose gas, in 
contrast to the Hubbard model of the fermion gas, the 
dilute limit describes already interesting physics such as 
BE condensation. The dilute limit of a Bose gas can 
be described quite adequately by traditional approaches 
such as the GP approach, which takes into account only 
the two particle interactions. However, it is intuitively 
obvious that as the density of particles in the BE con- 
densate increases, the nature of interparticle interactions 
will become more and more complex, and one will need to 
take even three-body and higher effects into account. It 
was demonstrated by one of us that the slave boson repre- 
sentation allows one to describe the dynamics of a Bose 
gas at arbitrary densities Q], which is what we review 
next. For the case of bosons, the slave boson represen- 
tation is even easier to formulate because there are only 
two states per site in a hard-core system: a site is either 
empty (represented by a complex field e x ) or occupied 
by a single boson (represented by a complex field b x ). A 
hopping process of a boson appears as an exchange of 
an empty site with a singly occupied site. Following the 
standard arguments |Q, this picture can be translated 
into a slave boson action S s .b. of the form 

X]E K e xtx.x'b x 'e*, + V x b* x b x + i\ x (e*e x + b* x b x - 1)]. 



Here again we have neglected fluctuations with respect to 
time because they have a gap fc^T at non-zero tempera- 
ture T as discussed in Sect. II. The field \ x enforces the 
constraint e*e x + b*b x = 1 which guarantees that a site 
is either empty or singly occupied. This becomes clear 
if we consider the partition function where we integrate 
over all fields 



Z= e- s °»Y[d\ x db x db x de x de x . 



(7) 



The A-field creates a Dirac delta function for the con- 
straint. The slave boson fields e x and b x can be combined 
to a collective field b*e x — ► & x . Then the constraint field 
A and the slave boson fields can be integrated out which 
leads finally to the action for the collective field <f>. This 
was demonstrated in detail in Ref. [Q. Here we only 
present the results: The new partition function reads 



Z = J e - s -- Sl Y[d<S> x d<S>* x 



(8) 



with hopping ("kinetic") term 

^> x {vH*) x 



6a 2 (l-r) 2 
and the potential term 



(9) 



-14/2-1/4 



df x e~ 



S inh(^/(<p x -V x /2) 2 + \<5> x \ 2 ) 
yJ{ Vx - V x /2f + |$,| 2 



(10) 



(6) 



Sb + Si is a generalization of the GP action to systems 
with arbitrary density. It agrees with Eq. (2) in the 
dilute regime after expanding Si in powers of \& x \ 2 up 
to second order: Si = ai|$ K | 2 + a^l^l 4 + o(|$| 5 ) can 
be compared with the potential of the dimensionless GP 
energy. For instance, we find for V x = the coefficients 
<zi = — 1/6 and 02 = 1/180. In general, both coefficients 
of the expansion depend on V x . That means the compos- 
ite field <3? of the slave boson approach must be rescaled 
by (813/2CI2) 1 / 4 in order to get the equivalent of the field 
$ of the GP approach. Furthermore, there is a renormal- 
ized chemical potential (1 — t) _1 — 1/6 instead of r in 
(2). t can be fixed by comparing the slave boson result 
and the GP result in the vicinity of a vanishing conden- 
sate. We found with the above mentioned parameters 
t w —5.5. Going beyond the dilute regime the effective 
potential of the slave boson representation deviates sig- 
nificantly from the GP potential: it grows only linearly 
for large in contrast to the |$ x | -behavior of the GP 
case. Therefore, the 'confinement' of the condensate is 
much weaker and allows a destruction of the latter by 
fluctuations. The strong potential in the GP case makes 
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the condensate very robust against fluctuations with suf- 
ficiently large \$ x \- 

An interesting quantity for the characterization of the 
condensate is the momentum distribution of condensate 
atoms H . For the case of a inhomogeneous condensate, 
the condensate wave function obviously is not an eigen- 
function of the momentum. However, we do expect the 
momentum distribution to be sharply peaked around zero 
momentum and its shape to be determined both by the 
trap parameters and the interparticle interactions . In 
the continuum limit that we are considering here, the 
momentum distribution can be evaluated from the of the 
field $ as § 

(*fc*lfc) = E eW " (a, "*' ) < $ -^>' C 11 ) 

where the average is taken with respect to the effective 
action (...) = Z^ 1 J ...e~ Sb ~ Sl J[ x d<£ x d$* x . In order to 
calculate this quantity we use a saddle point approxi- 
mation for the functional integral described in the next 
section. 



A. Mean-Field Theory 

The total density and the density of the condensate in 
the trap can be calculated again from the saddle point 
approximation. For the present case we feel that it is suf- 
ficient because thermal fluctuations will not be important 
at the small temperatures at which the BE condensation 
occurs in these systems. However, if one needed to ac- 
count for thermal fluctuations, one could do so by study- 
ing the deviations around the saddle point in the func- 
tional integral of Eq. (0) || . That means we have to look 
for a solution <& x which minimizes Sb + S\. This prob- 
lem is analogous to the minimization of the GP energy, 
discussed above. Instead of the non-linear Schrodinger 
equation Eq. (3), we have a generalization of this equa- 
tion for the slave boson approach 

ilO 

[(r/6a 2 )V 2 + (1 - r)]$ T + (1 - rf^ = 0. (12) 

This non-linear differential equation is even more compli- 
cated than the non-linear Schrodinger equation. There- 
fore, we apply again the Thomas-Fermi approximation 

r) 9 

^ + {l ~ T) d¥ x =Q - (13) 

dSi/d§* x contains higher order terms (three-body inter- 
actions etc.) which are important for the dense regime. 
However, inclusion of the kinetic energy term will only 
be necessary if more complicated structures (e.g., vor- 
tex states [[L0| ) are considered. Since the accuracy of the 
Thomas-Fcrmi approximation increases with increasing 
density, it is particularly suitable for our case where we 
are interested in the high-density regime. 



The total density of bosons can be evaluated from 
D\vlZ I 3[Lq. This quantity is measured with respect to 
a lattice Bose gas with lattice spacing a. The maximal 
density is n = 1, where we have one boson per site (i.e., 
one boson in a volume element a 3 ). In the trap the maxi- 
mum is at the center and it decreases monotonically with 
\x\ (cf. Fig. 1). 
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FIG. 1. Total density of bosons in a harmonic trap as a 
function of the distance \x\ from the center of the trap for 
the chemical potential fj,o = —0.6,0,0.6,1.6,2.6. The total 
density at a given distance |a;| increases as the chemical po- 
tential increases. This result assumes the equilibrium of the 
grand canonical ensemble of bosons, a condition which may 
be violated in the experiment. 

Solutions for an isotropic trap are presented for dif- 
ferent values of hq. In reality, of course, the number of 
bosons N is measured and /iq must be determined self- 
consistently for N while solving Eq.jT^). 
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FIG. 2. Density of the condensate in the trap, normalized 
by the volume, for the fio values of Fig.l. The surface of the 
condensate grows with increasing no, and the condensate for 
the lowest density has already the parabolic shape found in 
the GP approach. 

However, in this case, we have used /j,q as a parameter 
which controls the number N, and hence, the density of 
atoms in the trap and assigned it different values so as to 
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solve Eq.(|T^) for different values of particle density. So- 
lutions of Eq. ( |l3| ) for the density of the condensate are 
given in Fig. 2. The maximum of the condensate appears 
always at positions in the trap where the total density 
is 0.5. This reflects the duality of bosons and holes in 
the hard-core Bose gas on the lattice. If the total den- 
sity is less than 0.5 everywhere in the trap, the density 
of the condensate decays monotonically with \x\. This 
result demonstrates that the condensate is depleted at 
the center of the trap if the total density of bosons is 
larger than 0.5. Consequently, it is difficult to create an 
extended condensate in a potential where the bosons are 
concentrated at the center with high density. However, 
known experiments on a Bose gases are far away from 
such high densities. 

The correlation function of the order parameter field 
$x = l^zle 11 ^ can be approximated by neglecting the 
fluctuations of |<&J as 



(14) 



where <& x is a solution of Eq.(13). (The global phase of 
$ x , of course, is not determined by Eqs.(p^) or (|l^) due 
to symmetry.) The phase coherence of the fluctuating 
phase (p x may be larger than the size of the condensate 
because of the off-diagonal long range order in a three- 
dimensional Bose gas fl5|| . Therefore, the momentum 
distribution, defined by Eq . ( fil| ) and the continuum limit 
can be approximated by [pi 



/(*) 



d 3 xe lkx ^ a 



(15) 



The momentum distribution of the condensate atoms ob- 
tained using Eq. ( |i~5| ) , based on the solution of Eq. (|l^) is 
plotted in Fig. 3 for several values of the chemical poten- 
tial /io- 



IV. CONCLUSIONS 

In conclusion, we have studied the condensation of a 
three-dimensional high-density Bose gas in a harmonic 
trap. We demonstrated the unphysical nature of the so- 
lutions that one obtains if one applies the traditional ap- 
proach of Gross-Pitaevskii to study such a system. We 
identified that the reason behind this behavior was that 
the GP action ignores the three-body and the higher or- 
der interactions which become important at high densi- 
ties. In this work we proposed an alternative approach 
based on the slave boson representation, which accounts 
for these complex interactions at high densities by sat- 
isfying the hard-core condition strictly. Our approach 
leads to solutions which are well behaved at high densi- 
ties and predicts depletion of the condensate in the re- 
gions of high densities (mainly center of the trap) which 
one would expect on intuitive grounds. We also study 
the momentum distribution of the atoms in the trap and 
observe a narrowing of the momentum distribution in the 
high-density limit. This approach is expected to be of lit- 
tle application for the low-density condensates which are 
being created in the experiments presently, but we hope 
that it can be tested in future when experimentalists may 
be able to realize a dense Bose gas in a trap. We also 
plan to apply this approach to study the structure of vor- 
tices in a high-density inhomogeneous Bose gas, results 
of which will be presented in a future publication. 

Acknowledgement: We would like to thank Professor 
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FIG. 3. Momentum distribution of condensate particles in 
the trap for the yLto values of Fig.l. The sharpness of the 
distribution increases with increasing yuo 

Upon inspecting Fig. 3, one immediately observes that 
the momentum distribution becomes even sharper due to 
depletion. This is a consequence of spreading of the con- 
densate due to depletion which supports small momenta. 
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